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ABSTRACT 

The recent ASCA results show that the iron abundance of the X-ray gas 
of eUiptical galaxies is less than the solar abundance {Zq). The observed low 
iron abundance is inconsistent with the predictions of the previous chemical 
evolution models. In order to solve this problem, wc present a simple model of 
chemical evolution for elliptical galaxies after the galactic wind period under the 
assumption that the gases ejected from stars do not mix with the circumferential 
gas. The ejected gas components evolve separately according to their birth 
time and origin. We have investigated their evolution qualitatively. The gas 
components originated from supernova remnant shells cool and drop out of the 
hot gas faster than the other components because of their high density and metal 
abundance. As a result, supernovae cannot heat the whole gas of the elliptical 
galaxies effectively, in contrary to the previous results. If the metal abundance 
of mass-loss gas is not uniform, the mass-loss gas with higher abundance also 
easily drops out and the average abundance can decrease. We believe that this 
is a hint of solving the low abundance problem. 

Subject headings: galaxies : intergalactic medium - galaxies : interstellar matter 
- galaxies : X-rays 
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1. INTRODUCTION 

It is now well established by X-ray observations that elliptical galaxies are often 
sources of thermal X-ray emission. Temperatures of the hot interstellar medium (ISM) 
are well-determined to be around 1 keV and X-ray luminosities are typically around 10^° 
erg. In contrast, until recently, we have had little information on the metal abundance 
and the spatial structure of the hot ISM due to the poor energy and spatial resolution of 
detectors. For this reason, in most of works on the chemical evolution for elliptical galaxies, 
it had been assumed that the gas ejected from stars spontaneously mixes with the ISM. 
This inevitably leads to the metal abundance of the ISM higher than the value for the sun, 
because the metal abundance of the stars in elliptical galaxies is higher than the value for 
the sun (e.g. Arimoto & Yoshii 1987 ; Matteucci & Tornambe 1987 ; Mihara & Takahara 
1994 ; Matteucci & Gibson 1995 ; Fukumoto & Ikeuchi 1996). 

The recent observations of ASCA (Awaki et al. 1994 ; Matsushita et al. 1994) have 
changed the situation of this problem significantly. The ASCA observations showed that 
the iron abundance of the hot ISM is less than that of the sun, which is inconsistent with 
the predictions of the previous simple chemical evolution models. 

A hint of solving this problem is given by the recent ROSAT observation and 
theoretical argument. First, the recent i2(95'>lT' observation showed the inhomogeneity and 
the existence of cooling clumps in the hot ISM of NGC 507. (Kim & Fabbiano 1995) 
Second, Mathews (1990) indicated that the 1-^0 of metal ejected by each supernova 
event into the ISM is trapped locally within the hot bubbles and will never mix thoroughly 
with the ambient hot gas. He also argued that metal ejected through stellar wind may not 
mix either and make small fluctuation of metal abundance in the ISM. These observational 
and theoretical arguments suggest that the ISM is not well mixed, which is inconsistent with 
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the assumption of the previous chemical evolution models. If this is the case, we expect 
that gas components with higher metal abundance and density will cool faster than other 
components ; they go into cool phase and do not emit X-ray. This will reduce effectively 
the observed average abundance of the hot ISM. Renzini et al. (1993) has already explored 
the similar idea based on two-phase model and showed that the model does not reproduce 
the observed iron K line of NGC 1399, the cD galaxy in the Fornax cluster. Although 
the galaxies we investigated here are non-cD galaxies, we suggest that their two-phase 
model is too simple to explain the observation. In this paper, we investigate the chemical 
evolution of elliptical galaxies assuming that the gas released from stars does not mix with 
the ambient ISM ; we propose the multiphase model consisting of hundreds of gas phases 
and consider their evolutions. 

The paper is organized as follows. First, in the next section we describe the fundamental 
assumptions in detail and write down the basic equations to determine the evolution of 
the iron abundance of the ISM. Then, in §3, we solve these equations numerically for some 
typical models. On the basis of this general analysis, in §4, we simulate the X-ray spectra 
of our model galaxies and discuss the possibility of accounting for the observed low iron 
abundance of the hot ISM by comparing these simulations with the observational data. Our 
conclusions are summarized in §5. 

2. ASSUMPTION AND BASIC EQUATIONS 

We consider the chemical evolution after galactic wind period. We assume for simplicity 
that an elliptical galaxy is a sphere of radius R and it is uniform on average, although it 
contains many gas components (hereafter called "phases"). The phases are classified by 
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their origin and birth time. We consider the three types of the origin, that is, gas ejected 
through stellar wind (mass-loss gas), shell and cavity of supernova remnants (SNRs). From 
now on, the phases of each type are called "mass-loss phase" , "shell phase" , and "cavity 
phase" indicated by the indices ML, s, and c, respectively. We divide time into a finite 
number of steps, < to < < ••• < U < ••• < tn, where to and t„ are the start and end time 
of the calculation, respectively. We define t = to(= O.SGyr) as the time when the galactic 
wind stops. In each time-step, one phase is born for each type. The phase of type a {ML, s, 
or c) which is born in tj_i < t < tj is called the {i, Q;)-th phase and has its own temperature 
j'{i,a) ^ density p*^*'"\ metal abundance and mass M(*'°). Each ph ase radiates its 

thermal energy and evolves. We assume that energy transfer between the phases is worked 
only by pressure and ignore thermal conduction for simplicity. The phases are assumed to 
be in pressure equilibrium because sound crossing time of an elliptical galaxy is far shorter 
than its age (see Eq.(^)). Exceptionally, the phases whose temperatures become below 
7"crit(= lO^K) are not considered to be in the pressure equilibrium because the condition 
of the pressure equilibrium will broken down for the phases due to the high cooling rate, 
although it also depends on size of each gas blob composing the phases (Cowie, Fabian, & 
Nulsen 1980). They are assumed to cool immediately and drop out of the hot ISM. This 
means that while the total number of the type of the phase is unchanged, the number of 
the phase whose temperature is above T^rit can be different among the types at the given 
time. The gas left after galactic wind is also considered as a phase called "zero-phase" or 
(0, 0)-th phase. 



In § |2.1| , we describe the formation of the phases, and t = ti{l < i < n) unless otherwise 



mentioned. In ^.2|, we describe the evolution of the phases and the galaxy. 



2.1. FORMATION OF EACH PHASE 
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2.1.1. Mass-Loss Phases 



Here, we describe the formation of mass-loss phases. The relation between the mass of 
a star m(T), and its lifetime r, is given by 



logio M't) = 1-983 - 1.054^1ogioT + 2.52 , (1) 

where the mass is in units of solar mass (Mq) and the lifetime is in units of Gyr. 
(Larson 1974). Because we consider the chemical evolution after galactic wind stops 
{t > to = O.SGyr), this equation implies that we do not have to consider stars with mass 
larger than 2.85Mq in our calculation period. We assume that stars with mass in the range 
of 0.1 — 2.85Mq lose their masses by stellar winds ; we simply assume that the mass loss 
occurs instantaneously at the end of the life given by the fraction / of the initial mass of a 
star with m, 

for m < 0.7 , 

/ = \ 0.42m for 0.7 < m < 1.0 , (2) 

0.8 - 0.43/m for 1.0 < m < 2.85 , 
(Koppen & Arimoto 1991). Further, we assume that the temperature of the gas ejected by 
the stellar wind immediately becomes equal to the virial temperature of the galaxy T for 
simplicity. 

Under this assumption, the initial temperature and density of the mass-loss phase are 
then given by 

T(i'^L\ti) = f , (3) 

where // is the mean molecular weight (= 0.6), mn is the mass of the hydrogen atom, k-Q 
is Boltzmann constant, and P{t) is the average pressure of the ISM. The pressure P{t) is 



- 7- 



obtained by solving the evolution equations of the galaxy described in [|2]2. 



The gas and iron mass of the mass-loss phase at its birth time, M^^'^^^\ti) and 
Mp^^^\ti), respectively, are given by 



{i,ML) 
Fe 



{ti 



u 
u 



- (ML) 



{t)dt , 



(5) 



(6) 



where l'^^^ and Lp^^^ are the gas and iron mass ejection rates from stars, respectively. 
Since the time scale of star formation in an elliptical galaxy, which is typically 10^^^ yr, is 
short enough compared with the galaxy age (e.g. Arimoto & Yoshii 1987), we assume that 
the stellar system of the galaxy formed at t = simultaneously. Thus, Li^'^^'' and L^^c^^ 



given by 



dm{ 



T 



ML 



f 



dr 
dm{T 



(m(r = t))M,(0) 



(7) 



T=t 



dr 



0(m(r = t))M,(0) 



T=t 



where (pim), M^{t) and Zml are the initial mass function (IMF), the stellar mass of the 
galaxy and the iron abundance of the mass-loss gas, respectively ; the power of the IMF is 
taken to be 1.35 and the other variables are given by later. 

The density p(*'*^^)(ti), the temperature T(*'*^^)(ti), the mass M^''^'^\ti), and the iron 
mass Mp^^^^\ti) determined by above equations give the initial conditions of the evolution 
equations of the phases (see § p.2|) . 



2.1.2. Shell and Cavity Phases 



Next, we describe the formation of the shell and cavity phases. Since the progenitors of 
Type II supernovae have short lives t < to, Type II supernovae are important only in the 
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galactic wind period (e.g. Arimoto & Yoshii 1987 ; Matteucci & Tornambe 1987 ; David, 
Forman, & Jones 1991 ; Matteucci & Gibson 1995). Therefore, we consider only Type la 
supernovae (SN la) because we are concerned with chemical evolution after galactic wind 
stops. 

We use the adiabatic Sedov solution in pressure equilibrium to determine the initial 
conditions of the shell and cavity phases for simplicity. The actual shape of the SNRs is 
deviated from the Sedov solution because of the high temperature of the unshocked gas. 
However, we confirmed that the deviation is not important (see 

For simplicity, we assume that the SNRs are consists of two parts, outer shell and inner 
cavity region, and that they evolve as separate phases. We call the former "shell phase" and 
the latter "cavity phase", denoted by the indices s and c, respectively. When the distance 
from the center of the SNR is r and the shock front radius is Vg, the inner cavity region and 
the outer shell region correspond to r < (1 — k)rs and (1 — k)rs < r < Tg, respectively, where 
k is representing the width of the shell. In each region, the average density and temperature 
are decided as follows. 



The shock front radius is given by 

/9 09 Ro.x:^2 \ 1/5 

(9) 



VP/ 

(Spitzer 1978), where Esn{= lO^^erg) is the energy released from a supernova, ts is the time 
past since the supernova explosion, and p is the average density of the ISM. The pressure 
just inside the shock front is given by 

' (7 + 1) [dtj 25(7+1) ^ ' 

(Spitzer 1978), where 7(= 5/3) is the adiabatic constant. 

We assume that Pi is equal to the average pressure of the ISM P. Under the 
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assumption, the time tg is given by 

ts 

and 



-1 5/6 



25(7 + 1) 



1/3 



(11) 



(12) 



_25(7 + 1) 

The density and pressure distribution of the Sedov solution, psedov(^) and Psedov(^)) 
respectively, are determined uniquely under the condition P — Pi. 

In terms of and psedov(^)) the average density in the shell region is expressed by 

Ps = / ^'^Psedoviry dr/ / iirr dr . (13) 

J(l-k)ra J(l-k)rs 



The mass inside the shock front is given by 



where mpro is the mass of the progenitor star. The mass in the shell region is given by 



(14) 



rr 
J{1 



m^ — p^ I Anr dr . 

'(l-fe)rs 



(15) 



Hence, the mass and average density in the cavity region are given by 

rric = msNR - , 
3mr 



(16) 
(17) 



47r(l - k)^r^ ■ 

The iron abundance is assumed to be uniform inside the shock front radius Tg, because 
Mathews (1990) estimated that in an elliptical galaxy iron fragments made by Rayleigh - 
Taylor instability reach the shock front ~ 10^ yr after the supernova explosion. The iron 
mass in the cavity and shell region are 



VmsNR 



+ Z , 



(18) 
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mFe,c = rris 



i^ + z), (19) 



VmsNR 

where mpe is the iron mass ejected by the supernova and Z is the average iron abundance 
of the ISM occupied by the SNR. 

Next, we determine the temperatures. For the Sedov solution, the specific thermal 
energy in the shell region is given by 

U, = — r ^^Anr'dr . (20) 
nis J{l-k)r, 7-1 

The specific thermal energy in the cavity region is given by 

f/e = (^sN - mJJ,)/m, . (21) 

Since the thermal energy of the ISM occupied by the SNR cannot be ignored, the average 
temperatures of the shell and cavity region are given by 

r. = |^^!^ + t. (22) 

T. = \^ + f, (23) 

3 Kb 

respectively, where T(= Pfim-n/kBp) is the average temperature of the ISM which is 
assumed to be equals to the virial temperature. The densities and temperatures determined 
so far are adopted as the initial values of the shell and cavity phase for the evolution 



equations of the phases described in § |2.2| as follows : 

p'~'''\U) = p, , (24) 

T^'''\U) = Ts, (25) 

P^^'^H^.) = Pc , (26) 

r(^'^)(tO = T,. (27) 



- 11 - 



From Eqs. (|T2]) - (^), they are determined uniquely if -Esn, n^pro, f^Fe, k, p, T, and Z 
are given. We fix -Esn, f^pm, and mpe by giving typical values, that is, mpro = 2.0Mq, 
mpe = 0.75Mq, and i?sN = 10^^ erg. We also fix k = 0.3 ; the results in §^ are qualitatively 
unchanged even if we take the different value for 0.25 <^ k <^ 0.45. The other parameters, 
p, T, and Z are time-dependent and determined by solving the evolution equations of the 
galaxy described in § ^.2| . For examples, for p = 1.67 x 10^^''g cm~'^, T = 1.3 keV, and 
Z = I.OZq, we can estimate each quantity : ps = 2.4 x 10~^'^g cm~^, pc = 2.3 x lO^^^g cm~^, 
Ts = 1.8 keV, Tc = 9.6 keV, and mpe.s/'^s = ^Fe,c/^c = 2.8Zq. Note that since the 
cooling function is nearly independent of temperature and inversely proportional to metal 
abundance near T ~ 10''K (see Eq.(p7D), the cooling time is proportional to T/{pZ). Thus, 
the cooling time of the shell phase (oc Ts/{psiTLpe s/^s)) is generally shorter than that of the 
ambient medium (oc T/{pZ)). It causes selective cooling of the shell phase. This effect will 
be discussed in 

The gas and iron mass of the shell and cavity phase at their birth time are given by 

M^'''\U) = r Li'\t)dt, (28) 



M^'''\ti) = r Li'\t)dt, (29) 



i-l 



Mt\u)= r 4'ht)dt, (30) 

Mt\u)= r L^lm. (31) 

where Li^\ Lp^, and Lp^ are the mass and iron production rates of the shell and cavity 
phase, respectively ; they are given by 

LW(t) =rsN(t)m,(t), (32) 

Ll^\t) = rs^{t)m,{t) , (33) 

4i(^) ='^SN(t)mFe,s(t) , (34) 
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4e W ='"SN(t)mFe,c(t) , (35) 

where rsN(^) is the SN la rate. The time dependence of the SN la rate is assumed to be 
rsTssit) oc (David, Forman, & Jones 1990). The normahzation is given by the SN la 
rate at t = tn, which is specified later. 



2.2. EVOLUTION OF THE PHASES AND THE GALAXY 



The energy equation for the (z, a)-th phase for t > ti is given by 

/k,n.'>m\ _ kj>r^d ^ _(„..>)).A(2(.,.).r«.")w) . (36) 

7-1 dt\ uniH J /xmn dV ^ ' \ e ) \ \ )! ^ \ ) 

where n^*'") is the electron density of the (z, a)-th phase and A is the cooling function. 
Assuming that the relative metal abundance is the same as that of the sun for simplicity, 
the cooling function is approximated by 



2.1 X 10"^^ 1 + 0.1- 



1g 







+ 0.04 + 



6.2 X 10" 



-19 





(ergs cm~^s~"'^) 



'rp{i,a) - 

K 



-0.6' 



(37) 



This is an empirical formula derived by fitting to the cooling curves in Figure 9-9 of Binney 
& Tremaine (1987). 

In our model, all existing phases are assumed to be in pressure equilibrium because 
sound crossing time of an elliptical galaxy is far shorter than its age (see Eq. (^7]) ) . Thus, 
for i < j, 



Pitj. 



pmn 



p(^'°)(t,)^^5^^^ = P(t,) 



3/ ' 



(38) 
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M.)^^, (3.) 

where Mg and V are the total gas mass and volume of the galaxy, respectively. Here, we 
have approximated dP/dR by P/ R for simplicity. Note that V = (47r/3)i?'^. 

The total gas mass, the iron mass and the gas volume are the summation of those of 
each phase ; 

Mg(t)=E E M^''''\t), (40) 
MFe(t)=E E Mt\t), (41) 

v{t)=Y: E y^''''\t), (42) 

where 

As mentioned above, the phases whose temperatures are below Tcrit(= lO^K) are not 
included in the summation. 

Note that part of the SNRs is composed of pre-existing phases. As mentioned in 



§ |2.1.2| , we assumed that the SNRs evolve as new phases, then the masses of the pre-existing 
phases are reduced by the occupation by the SNRs. Thus, the mass and iron mass of the 
pre-existing (i, a)-th phase at t = tj [i < j) are given by 



(44) 



(45) 



t=t 



j-1 



-M(^'")(t) V{t) 

respectively, where Mg^^p^ is the mass occupied by SNRs during tj^i < t < tj, which is given 
by 

iV^^NR=r ^Mt)rs{tfrsMt)dt . (46) 
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We take the sound crossing time as the time-steps for the calculations 

3/2 / \ -1/2 



Riti) ..7 ( R V ( MiR) 

t - 1 . = ^ ^' = 4.7 X lOV ^ . 47 



Strictly speaking, calculations under Eq. ( ^Tj) are not exact because the phases whose cooling 
time is shorter than the time-step are not treated correctly ; these phases are removed at 
the next time-step tj^i and their PdV work in tj < t < tj^i is not taken into account. 
However, we confirmed that in our examined models their PdV work is less than 10% of 
the variation of the total thermal energy of the galaxy in tj < t < tj+i. Then, we believe 
that it is negligible for our qualitative discussion. 

We derive p^'-^H^i) and T^'''^\tj) from p(^'")(tj_i) and T^'''^\tj-i) by iterating Eqs.(|3|), 
(H), (H), (140), (E|), and until they converge. 



3. NUMERICAL ANALYSIS 



The main purpose of this section is to look into the qualitative features of the multiphase 
model of the hot ISM of the elliptical galaxy. We reduce the number of the free parameters 
by giving typical fixed values to some of them. First, for the parameters regarding the 
whole galaxy, we take to = 0.5 Gyr, = 10 Gyr, R{to) = 50 kpc, M*(0) = 2 x W^Mq, and 
Mg^Feito) / Mg{to) = 4Zq. We leave Mg(to) as a free parameter. Then if we specify Mg(to), 
we can write 1^(to)^°'°^ = (47^/3)i^(to)^ M(0'0)(to) = Mg(to), and M^l'°\to) = MFe,g(to)- We 
take the temperature of the zero-phase T^^'^\to) = 1.3 keV. The virial temperature of the 
galaxy is given by 

f = T(0'°)(to)(=1.3keV). (48) 
Here, we discuss the validity of adopting the Sedov solution as the initial conditions of 
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the shell and cavity phases by comparison with the numerical simulation of evolution of a 
SNR in an elliptical galaxy. Mathews (1990) calculated the evolution of a SNR expanding 
into ISM with temperature T = 10'' K and density h = 10~'^cm~'^ until the pressure 
equilibrium is achieved. Figures 3 and 4 in his paper depict temperature and pressure 
distributions when the pressure equilibrium is almost achieved [tg = 1.17 x 10^ yr). From 
now on, we call these the "Mathews solutions" . Although Mathews showed the variations 
of temperature and pressure with radius only in the case of T = 10^ K and n = 10~^cm~^, 
the fact that the radius where the gas temperature has risen by ~ 2 above the ambient 
value (T = 10'' K) varies as n^^^^ (his Table 5) may indicate that the solution is self-similar 
; we assume the self-similarity in the following argument. In the Mathews solution, the 
temperature and pressure gradients of the SNR are somewhat smaller than those of the 
Sedov solution. Thus, to obtain the similar values of ps, Pc, Ts, and Tc, the shell width k 
must be larger when the Mathews solution is used in Eqs. ([13D and ( ^OD than that when the 
Sedov solution is used. For example, when T = 10'' K and fi = 10~^cm~^, the values of ps, 
Pc, Tg and derived by using the Mathews solution and taking k = 0.35 are equal within 
30% errors to those derived by using the Sedov solution and taking k = 0.3. We confirmed 
that the variation of 30% does not affect the results in this section significantly. Moreover, 
we confirmed that the results in this section are not significantly changed even if we use the 
Mathews solution and take k = 0.3. 

We solve the basic equations described in §|] for the models whose details are given in 
Table 1. The difference between the models A and B is the initial gas mass Mg(to), which 
is chosen to match the observed X-ray luminosities of elliptical galaxies (~ lO^^ergs^^) at 
t = tn (Table 2). Models Al - A3 and Bl - B3 are calculated to see the effect of SN la. 
The SN la rate is normalized by the values at t = t„(= lOGyr) shown in the third column 
of Table 1. The values are expressed in units of SNu, that is, the number of SNe per 
lO^^/i^^LgQ per 100 yr [Hq = lOO/ikms^^ Mpc^^ ; we set h=0.5). The ratio of the mass to 
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luminosity is taken to be 8Mq/ Lq, and the rate for SN la is normalized by using this ratio. 
Models A4 - A5 and B4 - B5 are calculated to see the effect of metal- abundance distribution 
for the mass-loss gas. In these models, the mass-loss phase born in each time-step is 
divided into two phases with equal masses and different abundances corresponding to the 
two figures for Zml (Table 1). On the other hand, in models Al - A3 and Bl - B3, each 
component has the same abundance or the mass-loss gas has one kind of the abundance. 
Table 1 also shows total number of the time-step n and the numbers of phases which have 
survived (their temperature is above Tcrit) until t = tn- In models Bl - B3, all phases cool 
and vanish by sometime before t = tn, and we stopped the calculations at that time. We 
give the numbers of the survived phases in columns A^ml.h, A^ml.l, A'sh, and A^cav denoting 
respectively the mass-loss phase with the high and low abundance, shell phase, and cavity 
phase (Table 1). Figures 1-4 show the time evolution of R, Mg, and iron abundance Z, 
where Z represents Mg^pe/Mg in units of the solar abundance Zq, 1.7 x 10^^. 

Figures 1 and 2 show the influence of SN la on the gas evolution. The sharp decreases 
at t ~ 4 Gyr in Figure 1 and at t ~ 3 Gyr in Figure 2 represent cooling and vanishing of 
the zero-phase. Since the initial volume in series A is the same as in series B, the initial gas 
density is higher in series B than that in series A. Therefore, the cooling process is more 
effective and R decreases faster in series B than in series A. The radius of the gas sphere 
does not significantly increase even if the SN la rate at t = t„ increases (Fig. la and 2a). 
Furthermore, the gas mass decrease faster in the models of higher SN la rate (Fig. lb and 
2b). These facts indicate that supernovae are not effective heating sources of the ISM as a 
whole. The reason is that the shell phases cool and drop out faster than other phases owing 
to their high density and metal abundance (see § p.l.2 ), and that they have large thermal 



energies when they are born as the following arguments shows. When a SNR becomes in 
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pressure equilibrium, the thermal energy within a shell region is given by 



rUsTs + rricTc 



rUsT, 




2TikT 




(49) 



Since Eg ^ E'sn in our parameter range, more energy than that supplied by supernova 
explosion is radiated from the shell region in a short time. Thus, supernovae cannot heat 
the ISM. This is contrary to the previous results (e.g. Ciotti et al. 1991 ; David et al. 
1991) which claimed that supernovae are effective as heating sources and swell the ISM. 
However, this fact does not reject the galactic wind model. While the expanding SNRs do 
not overlap in our calculations, it is expected that they overlap during the galaxy-formation 
period because of the high supernova rate ; then, the ISM is expelled and galactic wind 
occurs (Ikeuchi 1977). Note that Mathews (1990) indicated the possibihty of buoyant 
mixing of the hot cavity region in the local ISM. This may decrease entropy fluctuation and 
the effect of selective cooling of the shell phases. 

In models Al, A2, Bl, and B2, Z is around Zq for i ;^ 4 Gyr (Fig.lc and 2c), because 
most of the mass of the remaining gas is dominated by mass-loss phases {Z = IZq). In 
models A3 and B3, the shell phases are dominant in mass ; they start to vanish for t ^6 in 
model A3, and for i 4 in model B3. 

As seen in Figures 3(c) and 4(c), Z is small when the mass-loss gas has two kinds of 
the iron abundance. This is because the mass-loss phases with the higher iron abundance 
also drop out in addition to the shell phases, as shown in Table 1. In models A4, A5, B4, 
and B5, Z < Zq at t = 10^° yr, which is consistent with the observations (Awaki et al. 
1994 ; Matsushita et al. 1994). Although our model may be too simple to give quantitative 
predictions, we believe this is a hint of solving the low abundance problem of elliptical 
galaxies. 
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4. SPECTRAL SIMULATIONS 



In most of observations, X-ray spectra of galaxies are fitted to at most two thin thermal 
plasma models. However, if the actual ISM is composed of many components with various 
temperatures and metal abundances, this simple fitting may give deferent values from the 
actual ones. In this section, we examine this possibility by constructing composite spectra 
from the many phases which have survived at t — tn{— 10 Gyr) and fit them to a few 
spectral models by the XSPEC package (ver. 8.50) as most of observers did, although our 
model may be too simple for detailed discussion. 

We examine models Al, A2, and B4 because in the other models X-ray fiux is below 
the observational limit (Table 2). The distance to the model galaxies is assumed to be 15 
Mpc. In our models, several hundred phases are survived until the present epoch. However, 
the XSPEC package cannot treat so many phases ; for each type we collect each 10 phases 
into a new single phase. We regard the luminosity-weighted temperature and abundance of 
the 10 phases as the temperature and abundance of the new single phase. We confirmed 
that even if we collect each 20 phases into a new single phase, the results do not change 
significantly. First, the spectrum of each new phase is calculated by a thin thermal plasma 
model (so-called Meka model ; Mewe, Gronenschild, and van den Oord 1985 ; Mewe, Lemen, 
and van den Oord 1986 ; Kaastra 1992) in the XSPEC package. The reason of using Meka 
model is that the XSPEC package can treat high- abundance components only with this 
model. Then, using the XSPEC package, the spectral components are composed and the 
Galactic absorption (lO^^cm"^) is taken into account. Since the cosmic iron abundance to 
which we referred as the solar abundance is 2.6 x 10~^ in the XSPEC package, we corrected 
it to 1.7 X 10~^ at composition. The model spectra are fitted to the Meka model plus an 
absorption with a response function of ASCA SIS. Background is not considered. 
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The exposure time and the results are shown in Table 2 ; the fluxes do not include 
absorption. For models Al, A2, and B4, the upper lines are the results of the fitting to 
one Meka model and the lower lines are those to two Meka models. For comparison, Z at 
^ = ^n(= 10 Gyr) is also shown. The abundance Zml, temperature Tst, and absorption 
column density A^H,fit derived by the fitting are also shown in Table 2. As clearly seen 
in Table 2, Zfit is much smaller than Z when the model spectra are fitted to one Meka 
model plus absorption. However is so large that these models will be rejected ; the fits 
suffer from significant positive residuals above ~ 2 keV as exemphfied in Figure 5. We 
next fitted the spectra with a sum of two Meka models with different temperatures and 
normalizations but using common abundance as Matsushita et al. (1994) did. In these 
cases, is significantly reduced. The fitted spectra are shown in Figure 6. All the models 
examined here are acceptable and the low abundance, Z^^ < Zq, is consistent with the 
observations (Awaki et al. 1994 ; Matsushita et al. 1994). The derived abundance Zfn 
is still smaller than Z. This may indicate that not only statics but also estimation of 
equivalent width of the lines is required to derive the metal abundance of multiphase gas 
(Matsushita, private communication). The iron K line is very weak (Figures 5 and 6), 
because of the small number of the survived shell phases (Table 2), and of the low density 
and high temperature (~ 10 keV) of the cavity phases. The weak iron K line is consistent 
with the observations (Awaki et al. 1994 ; Matsushita et al. 1994). 

5. SUMMARY AND CONCLUSIONS 

In this paper, we have presented a model of the chemical evolution for elliptical galaxies 
after galactic wind period under the assumption that the gases ejected from stars do not 
mix with the circumferential gas. The ejected gases evolve separately according to their 



-20- 



birth time and origin. We considered three origin of the ejected gas, that is, shell and cavity 
of supernova remnants and mass-loss. We fitted the model spectra made from the survived 
phases by thin thermal plasma models. 

The main results and conclusions can be summarized as follows: 

(i) The model predicts that the supernovae are not effective as heating sources of the 
ISM after the galactic wind stops. The shells of the supernova remnant can cool and drop 
out rapidly because of their high density and metal abundance. Since the shells initially 
have large thermal energy, the energy ejected by supernova explosions is radiated and does 
not heat up the ISM. 

(ii) The present iron abundance of the hot ISM can be less than IZq and consistent 
with the ASCA observations if the abundance of the mass-loss gas has wide distribution. In 
this case, the mass-loss gas with higher abundance also drops out in addition to the shells 
of the supernova remnant. This effect makes the observed average abundance of the ISM 
lower. Although our model may be too simple to give quantitative predictions, we believe 
this is a hint of solving the low abundance problem of elliptical galaxies. 

(iii) If the model spectra are fitted by one thermal spectrum model, the resultant 
abundance is well below IZq, although the fits suffer from significant positive residuals 
above ~ 2 keV. When they are fitted by two thermal spectrum models, is significantly 
reduced and the derived abundances are also subsolar. 
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Figure Captions 



Fig. 1. — The evolution of (a) radius (b) gas mass and (c) abundance for models Al - A3. 

Fig. 2. — The evolution of (a) radius (b) gas mass and (c) abundance for models Bl - B3. 

Fig. 3. — The evolution of (a) radius (b) gas mass and (c) abundance for models Al, A4 
and A5. 

Fig. 4. — The evolution of (a) radius (b) gas mass and (c) abundance for models Bl, B4 and 
B5. 

Fig. 5. — Simulated X-ray spectra observed with the ASCA SIS. The hne shows the best 
fitting (one Meka plasma + absorption column density), (a) Model Al (b) Model A2 (c) 
Model B4 

Fig. 6. — Simulated X-ray spectra observed with the ASCA SIS. The line shows the best 
fitting (two Meka plasma -|- absorption column density), (a) Model Al (b) Model A2 (c) 
Model B4 
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TABLE 1 
Parameters of Models 



Model 


(lOi"M0) 


SN rate 
(SNu) 


(Zq) 


n 


A^'ml.h 


A^ML,L 




^cav 


Al 


0.65 


0.05 


1.0 


133 


91 


91 


10 


133 


A2 


0.65 


0.02 


1.0 


138 


33 


33 


7 


138 


A3 


0.65 


0.2 


1.0 


117 


107 


107 


23 


117 


A4 


0.65 


0.05 


1.5,0.5 


128 


31 


126 


25 


128 


A5 


0.65 


0.05 


1.8,0.2 


124 


37 


124 


33 


124 


Bl 


0.9 


0.05 


1.0 


147 










B2 


0.9 


0.02 


1.0 


136 










B3 


0.9 


0.2 


1.0 


161 










B4 


0.9 


0.05 


1.5,0.5 


213 


21 


49 


12 


213 


B5 


0.9 


0.05 


1.8,0.2 


168 


24 


153 


21 


168 



1 



TABLE 2 
Results of Fitting 



Model 




Flux 


Exposure 


Z 








xVdof 




(lO-^Oergs-i) 


(10~-^^erg cm~2 s~-^) 


(ks) 


(%) 


(%) 


(keV) 


(102icm-2) 




Al 


24.6 


9.1 


40 


1.1 


0.03 


0.25 


0.26 


200/35 












0.65 


0.23, 1.0 


0.08 


32.4/32 


A2 


13.6 


5.1 


40 


1.0 


0.10 


0.32 


0.11 


48.1/27 












0.36 


0.33, 1.2 


0.00 


32.4/24 


A3 


1.5 


0.5 




2.3 










A4 


0.3 


0.1 




0.8 










A5 


0.1 


0.03 




0.7 










Bl 


















B2 


















B3 


















B4 


8.6 


3.2 


80 


0.8 


0.04 


0.41 


0.20 


133/42 












0.43 


0.40, 1.0 


0.02 


37.5/39 


B5 


1.5 


0.6 




0.5 











2 



